function x = UnitLxbSol(L,b,n)
%UnitLxbSol 解单位下三角阵方程
%输出结果x储存在b中
for j = 1:n-1
    b(j+1:n) = b(j+1:n) - L(j+1:n,j)*b(j); 
end
x = b;
end

